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The general treatment of a separable Hamiltonian of Liouville-type is well-known in operator 
formalism. A path integral counterpart is formulated if one starts with the Jacobi's principle of 
least action, and a path integral evaluation of the Green's function for the hydrogen atom by Duru 
and Kleinert is recognized as a special case. The Jacobi's principle of least action for given energy 
is reparametrization invariant, and the separation of variables in operator formalism corresponds to 
a choice of gauge in path integral. The Green's function is shown to be gauge independent, if the 
operator ordering is properly taken into account. These properties are illustrated by evaluating an 
exact path integral of the Green's function for the hydrogen atom in parabolic coordinates. 
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■ I. INTRODUCTION 

The path integral treatment of the hydrogen atom is interesting not only for a methodological 
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interest but also for a pedagogical purpose. After several early attempts [1,2], Duru and Kleinert[3] 
showed an elegant path integral method to evaluate the Green's function for the hydrogen atom 
exactly. Two basic ingredients in their method are the use of a re-scaled time variable and the 
| so called Kustaanheimo-Stiefel transformation [4] which reveals the 0(4) symmetry explicitly in the 
coordinate space. Many of the clarifying works of this approach have been published [5] - [14]. 
^ ■ It has been shown elsewhere [15] that this problem is treated in a more general setting of gauge 
theory if one formulates the problem on the basis of the Jacobi's principle of least action, which is 
reparametrization invariant. 

The general technique of gauge theory is thus applicable to the evaluation of path integral, and a 
suitable choice of gauge simplifies the problem such as the hydrogen atom. In particular, the Green's 
function is shown to be gauge independent. The use of the Kustaanheimo-Stiefel transformation is 
rather technical and it is not essential in solving the problem exactly. We in fact show a simple trick 
in parabolic coordinates [16] which solves the hydrogen atom exactly. 

. ^ ' II. HYDROGEN ATOM 

X: 

2.1, Analysis in Parabolic Coordinates 

We analyze the hydrogen atom by starting with the Hamiltonian written in terms of parabolic 
coordinates 

where the parabolic coordinates are introduced by £ = |(r — z), r) — |(r + z) and ip stands for the 
azimuthal angle around the z axis. We further perform a canonical transformation which simplifies 
the kinetic term in H as 
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M 2 , < U < OO 



77 = -f 2 , < v < oo 



and the Hamiltonian becomes 



2m u + v z u l v v z v u l + f 2 

where r = £ + r\ = {u 2 + t> 2 )/4. This Hamiltonian is not yet a separable one of Liouville-type. 
One may solve the Schroedinger equation 

Hi) = Etp 

or equivalently 

H T i) = 

with 

Ht = 7^lPl + ^Jl +Pl + ^V%\ -e 2 + ^f(u 2 + v 2 ) 
where ui is defined by 

-mu = — E 
2 4 

We consider the case E < for the moment. Ht stands for the total Hamiltonian defined by 
specific gauge condition; a general definition of Ht will be given later in (39). 
Eq.(5) may be rewritten in an equivalent form as 

H T ^J = 

We here introduced auxiliary variables (£y,<//) as 

Ht = ^[fu + ^Jl +PI + + ^(u 2 + v 2 ) - e 2 

= —f u + —u 2 + ±-? v + - e 2 

2tt/ u 2 2771 * 2 



and we defined 



■u = (m 1 ,-u 2 ) = (it cos ip, u sin y?) 

-2 -2,^-2 
Pu=Pu+ —P V 

v — (vi,V2) = (v cos ip' , v sirup') 

-0, _ -2 i 1 -2 /- 
Pv — Pv ' ^2^V V - 



The subsidiary condition in (8) replaces the use of the Kustaanheimo-Stiefel transformation, and at 
the same time it renders a Hamiltonian of Liouville-type. This introduction of auxiliary variables 
(10) has been discussed by Ravndal and Toyoda[16]. 

A general procedure to deal with a completely separated operator such as Ht in (9) is to consider 
an evolution operator for a parameter r defined by 

(u b ,v b \e- tilTT/h \u a ,v a ) 

= e^(u b \exp[-(i/h)(±? u + ^fu 2 )r]\u a ) 



1 _o muj 



x(v b \exp[-(i/h)(—f v + 



2m rv 2 



V )T\\Vn 



2irih sin ujt 

xexp{ ^ mUJ [{u 2 h + v b +u 2 a + v 2 a ) cos ujt - 2u b u a - 2v b v a }} (11) 
Aft sin ujt 

where we used the exact result for a simple harmonic oscillator [17]. 
A crucial observation here is that p^ and are preserved during the evolution dictated by the 

operator H T in (9), since [p^,, H T ] = [P<p',H T ] = 0. It is then sufficient to impose the constraint (8) 
only on the initial state , for example. Starting with a general state belonging to the eigenvalues 
p^ — m and = m' 

e im V mV (12) 

we can use the following trick 

Jo 2tt 

to project out the state satisfying p^ = p^r, and ip + ip' is regarded as the actual azimuthal angle. 
We thus obtain 

(u b , v b , (<p + if') b \e~ iHTT/h \u a , v a , (ip + <p') a ) 

= e (tt^T- / 7T ex P{^r- K + % + % + O cos^r - 2u b u a - 2v b v a }} 

2mksmujT Jo 2tt 2ft sm ujt 



^2nihsmujT Jo 2tt 2ft sin ujt 

e ( ? ) x 

2irift sin ujt 

r 2n d6 , imuj 



I 7T ex P{ [ 4 (£« + £b + Va + Vb) cos ujt - 4V2(r a r b + x a x b ) 1/2 cos(6» + 7)]} 

Jo 2n 2ft sm ujt 



= e 



ie 2 T . muj 2 2imuj , . 2\f2imuj _ _ 1/2 

L . fc ■ ) exp{— (r a + r b ) cos ujt \I (—. (r a r b + x a x b ) 1 ) 

2mft sin ujt ft sin cur ft sin cur 

(14) 

In this evaluation we start with the relation 

u b u a + v b v a = u b u a cos A<p + v b v a cos Aip' (15) 
with Ay? = ip b - ip a , Ay?' = (fi' b - (fi' a , and 
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UbU a cos(A<£> + 6) + vi,v a cos(A<// — 9 
= (ubU a cos Aip + f bV a cos A<//) cos 6* 
+(— Uf,u a sin A<^ + V(,v a sin Ay/) sin 6* 



= 4^ + Wa + 2\/tbLVbVa cos(Av? + Ay?') cos(# + 7) 
= 2\/2^/ r a r b + 2^ + p a p b cos(A<^ + Aip') cos(6 + 7) 

= 2\/2\J r a r b + x a x b cos(6» + 7) (16) 

where we used the definition of variables in (2), and 7 is a number independent of 6. We also defined 
a modified Bessel function 

2V2imu 1/2 [ 2 * d6 2s/2imu)(r a r b + x a x b ) 1/2 

Io{-r—. V t + x a x b ) 1 ) = / — exp{ — cos6>} (17) 

KsmuT Jo 2n ftsinwr 

The parameter r is arbitrary, and we eliminate r to obtain a physically meaningful quantity by 

i / dr(u b , v b , (ip + ip') b \e~ lHTT/n \u a , v a , (<p + </?%) 
Jo 

= (u b , V b , {if + <ff)b\-*r-\u a , V a , (<£> + 

. f°° , ie 2 T , mw 2 r 2?mw 

= W are ( - , ) exp{— {r a + r b ) cosur} 

Jo 2-KihsmuT asm our 

T ,2\f2imijj _ _ 1/2 

a sm cur 

mcu f 1 ,..^ 1 f -2muj 1 + A., r Ay/2mw A 1/2 _ _ 1/2 

= io rfAA (T^AF e ^ [ ^ (ra + rfc)( rrA )]/o( ^^T3A (rarfc + XaXb) } 

(18) 

where we rotated r by 90 degrees , r — > — rr, and defined 

A = e- 2u)T 

v = e 2 /2uj (19) 

We next show that (18) gives an exact Green's function for the hydrogen atom by noting the 
sequence 

(u b ,v b ,(pb\-t-\u a ,v a ,(p a ) 



= (Zb,Vb,<Pb\ J— ~ -\L,Va,(Pa)(—^ ) 

:( 77- , (£b,Vb,<Pb\£a,Va,<Pa) 



r]b, <p b ) - E V6 + Vb + Va 

471 p /2m - e 2 /r - E 



-^-(Xb\- 2 K) (20) 



where we used tp in place of tp + tp' and the relation (AB) 1 = B 1 A l . 
The volume element changes in this transition from to H as 

dVo = 2iruvdudvdip 
-> dV = (f + 77)dV{, = 4tt x 2(f + r])d^dr]dtp 

= Att x r 2 drd cos 9 dtp (21) 

The bra- and ket- vectors in (20) are normalized in the combination 

J dVo\u, v , ip) (u, v, tp\ = 1 
/ dV\£,ri,<p)—^—(£,r],<p\ = 1 

J t + T] 

d 3 x\x)(x\ = l (22) 



and the extra factor of An in dV = Anr 2 drd cos 6 dip requires the appearance of the factor of 1/An 
in the last expression in (20). The appearance of 2ir in dV is an artifact of the variable <p' in (10). 
This normalization condition of bra- and ket- vectors together with a symmetry in x a and Xb justify 
the identification (20). A more explicit and concrete analysis of eqs.(20)~ (22) will be given in 
connection with the Jacobi's principle later. 
As for the operator ordering, the momentum operator changes in (20) as 

Pi +fv = (^) 2 [^d u ud u + ^d v vd v ] 
= (-) 2 [^ + ^] 

= P&Pt + PvVPv ( 23 ) 

and 

1 1 ^2 

(-j-^)(fctfc + PvVPv) + ~^-P% = P ( 24 ) 

where the right-hand side is written in cartesian coordinates. We note that dV and dV in (21) 
respectively render Ht and H(£,r),ip) hermitian. 

Combining (18), (20) and (24), we have thus established the exact Green's function including the 
operator ordering 

(£&It9 \x a ) = ^r£- I d\\~ v 1 exp[ 2 ™ W (r a + r 6 )(j-^)] 

p/2m-e 2 /r - E nh Jo (1 — A) 2 h 1-A 

x/ (^ — Y^x^ aTb + a b) ] ( ] 

It is known that this formula, which was first derived by Duru and Kleinert[3], is a Fourier transform 
of Schwinger's momentum space representation[18]. The continuation to the scattering problem with 
E > is performed by the replacement 

uj — > (—i)cu, v — > iv (26) 

in the above formula. 
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One can understand the spectrum of the hydrogen atom by looking at H T in (9) [16]. If one defines 
the oscillator variables 

1 I moo i 
1 r fmw i 



one obtains 



ik = W— Vk+ v^i K] ' k = 1 ' 2 (27) 
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H T = hu[^2(a* k a k + a\a k ) + 2] - e 2 

k=l 

p v = ih[a\a2 — a 2 a i] 

p v > = ih[a\a 2 — a\ai] (28) 



After a unitary transformation 



ai = -^(h - ib 2 ) 

a 2 = ^(-ib 1 + b 2 ) (29) 
and a similar transformation of a\ and a 2 , one obtains 

H T = huj[j2(blb k + blb k ) + 2]-e 2 

k=l 

p v = h[b{b! - b\b 2 ] 

p v , = hlblh - b\b 2 ] (30) 
By defining the number operators 

n k = b{b k , 

n k = ~blb k , k = l,2 (31) 

the total Hamiltonian is given by 

Ht = hu[ni + n 2 + hi + h 2 + 2] — e 2 

= huo[2ni - + 2h 2 + p^/% + 2] - e 2 

= 2ftcj[ni + h 2 + 1] - e 2 (32) 

by noting the physical state condition p v = p^i. 
We thus define the principal quantum number (or its operator) n by 

n = n x +h 2 + 1 = 1,2,3, (33) 

and the physical state condition 

(2nhu - e 2 )ipp h y S = (34) 

gives rise to the Bohr spectrum 
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E = -^f(f) 2 - 2 , n = 1,2,3,... (35) 
2 he n z 

by noting the definition of uo in (7). As for the degeneracy of states with a fixed n, it is confirmed 
that each state with a given n has n 2 degeneracy. 
2.2, Jacobi's Principle of Least Action 

The meaning of the total Hamiltonian in (6) becomes transparent , if one starts with a Nambu- 
Goto-type Lagrangian ( the Jacobi's principle of least action for a given E ) which is reparametrization 
invariant, 



S= f T Ldr = f T dT\l2m(E-V(r))A 2 
Jo Jo V dr 



>2m(E -V(r))(dx) 2 (36) 
The momenta conjugate to coordinates are then defined by 



'"=§ = V»»<*-V(r))(£)/^)» (37) 

and one obtains a vanishing Hamiltonian as a result of reparametrization invariance and a first-class 
constraint <fi, which is the generator of reparametrization gauge symmetry, 

H =pS-L = Q 

-2 

<f>(x,p) = — + V(r) - E ~ (38) 
2m 

Following Dirac[19], the total Hamiltonian is defined by 

H T = H + a(x,p)(f)(x,p) 

= a(x,p)(j)(x,p) ~ (39) 

and the function a(x, p) specifies a choice of gauge and fixes the arbitrary parameter r in (36), which 
parametrizes the orbit for a given E. A change of the parameter r to r — 5/3(t, x,p) is generated by 
6/3(t, x,p)a(x,p)(f)(x,p) = 8(3(t,x,p)H t , for example, 

Sx(t) = x{t) — x(t) 

= x(t + 513) - x(t) 
= {x,5(3H T } PB 

= 5(3(T,x,p)^-x(r) (40) 

in terms of the Poisson bracket, since x'(j — 5(3) = x(t). 
Quantization is performed by 

ih^-4> = Ht^j (41) 
or 

with a physical state condition 

a(x,p)4>(x,p)i) phy = (42) 
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A specific choice of the gauge a(x, p) = r = £ + 77 leads to the Hamiltonian H T in (6) and the choice 
a(x,p) = 1 gives the original static Schroedinger equation(4), since the states ifi in (5) and (4) are 
physical states. Eq.(41) gives rise to the evolution operator in (11). 

We now explain the relations (20)~(22) in a more concrete manner. We start with eq.(20) for a 
generic negative E 

G(E; £ 6 , T] b , (p b ; £ , r] a , (p a ) = (£ b ,Vb,Vb\ ~ —77 -\£a,Va,<Pa) 

H T (E;£,r],ip) 

= J2 £ b , rj b , ip b )——(f)* n (E; rj a , <p a ) (43) 

with 

H T (E; £, rj, <p)(j)n(E; £, rj, tp) = \ n (E)<j> n (E; £, rj, tp) 
J <j>* n {E- £, 77, ip)ME; £, V, <P)dV = S n , t 

\ n (E) = 2nh ] J^--e 2 

dV = 4tt x 2d^drjdip (44) 

where we used the result in (32) and also the variables (£, 77, tp) instead of (u, v, tp) for notational 
simplicity. The summation over n in (43) is formal including the n 2 degeneracy. Note that the 
complete orthonormal states {<p n } in (43) are all unphysical off- shell states. In path integral, the 
summation in (43) is exactly evaluated in (18). 

We next rewrite G(E; rj b , ip b ; £ a , ?7 a , ip a ) in terms of physical on-shell states by writing an unsub- 
tracted dispersion relation (i.e., paying attention only to the pole structure in E) as 

G(E;£ b ,ri b ,ip b ;£ a ,Tl a ,ip a ) =J2M E n]^b,Vb,^b)— ^ 8\ n (E n ) . 4>n( E n,€a,Va,<Pa) (45) 

n {E n -E){ qe^~) 

by noting 

X n (E) = X n (E n ) + (E- E n )^^- = (E- E n )(^-) (46) 

for E ~ E n . 
When one defines 

i> n {E n ] £, 77, tp) = , == (f) n (E n ; £, r], <p) (47) 
V §K~ 

one can show the orthonormality relations of physical on-shell states 

J r n (E n ,^ V ^)MEi;^V^)^ + v)dV = J r n (E n ,^v^)MEi^,V,^)dV 

= 5 n ,i (48) 
with dV — (£ + r])dV . First of all, from the physical state condition 

H T (E n ;£,ri,(p)<f> n (E n ;Z,ri,(p) (49) 

1 111 

= {^{P&h + PvVPv) + g^(| + - ^n(£ + ?7) - e 2 }0„(£„; £, 77, y>) = 
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one can establish the orthogonality relation 



(E n -E t ) J <K(E n ^,r,,(p)MEi^,V,<P)^ + v)dV o = (50) 

for n 7^ /. Also from the relation (49) and the fact that the "eigenvalue" e 2 is equally distributed for 
the kinetic and potential terms for harmonic oscillators (in terms of u and v variables) , we have 

2 

-E n J ^ n (E n , £, 77, v)ME n , £, V, <f) (f + vWo = C - (51) 

namely 

dt(E n ) [ <t>n(E n \Z,y,<p)ME n \Z,V,v)(Z + v)dV = 1 (52) 

dE n 

by noting (44). This proves (48). 
From (45), we finally arrive at the expression 

h 

G(E] T) b , ^b] fa, Va, <Pa) = ^n{ E n\ 6>, Vb, Vb)^ ^ni E n\ Ca, Va, <Pa) 

n &n — & 

= (Zb,Vb,(Pb\-K— " -\£a,Va,¥a) 

H(t,V,<P) ~ E 
1 Ti 

= T~ (^Ito l^o) (53) 

47r p /(2m)-e 2 /r-E 

which establishes the gauge independence of the Green's function for negative E. Although we here 
used the same notation for the state |£, 77, V 7 ) in (43) and (53), the meaning of these states are quite 
different. This difference is explicitly exhibited in (20) ~ (22). It is important to realize that the 
exact path integral is performed for the off-shell states in (43). 



III. CONCLUSION 



The path integral treatment of a general separable Hamiltonian of Liouville-type is formulated on 
the basis of the Jacobi's principle of least action by using a gauge theoretical technique, and it has 
been illustrated for the hydrogen atom in parabolic coordinates. 

I thank C. Bernido, P. Cvitanovic, J. Klauder, L. D. Faddeev, V. Sayakanit, and H. Sugawara for 
helpful comments. 
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